Spontaneous emission in the presence of a spherical plasmonic cloak 
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We investigate the spontaneous emission of a two-level atom placed in the vicinities of a plasmonic cloak 
composed of a coated sphere. In the dipole approximation, we show that the spontaneous emission rate can be 
reduced to its vacuum value provided the atomic emission frequency lies within the plasmonic cloak frequency 
operation range. Considering the current status of plasmonic cloaking devices, this condition may be fulfilled 
for many atomic species so that we argue that atoms with a sufficiently strong transition can be used as quantum, 
local probes for the efficiency of plasmonic cloaks. 



I. INTRODUCTION 



Among the most interesting effects of Cavity Quantum 
Electrodynamics (Cavity QED) are the possible influences 
that bodies in the vicinities of atomic systems may exert in 
their radiative properties. The fact that the environment can 
strongly modify the spontaneous emission (SE) rate of a sys- 
tem was theoretically predicted for the first time (to the best 
of our knowledge) by Purcell in 1947, in a national confer- 
ence of the American Physical Society |[l|]. Approximately 
two decades later, Drexhage et al 10] observed the variation of 
the fluorescence decay time of a molecule near a mirror Also, 
Morawitz |3f] obtained the SE rate of a two-level system near 
a perfectly conducting mirror and showed that the result ex- 
hibited oscillations with the distance between the system and 
the mirror, showing that the SE rate can be enhanced or sup- 
pressed. The first experiments that confirmed these spacial 
oscillations for were done again by Drexhage in 1970 ]|3] and 
in 1974 ]jj]. Generalizations for an atom between two parallel 
conducting plates were done by many authors ll6l- [Tlll . and the 
possibility of modifying the emission properties of atoms and 
molecules in a controlled way has stimulated the investigation 
of SE in different systems and geometries [fl2ljl6ll . Even in 
simple situations as those with an atom between two parallel 
mirrors, quite interesting phenomena may occur. For instance, 
if the transition dipole moment is parallel to the mirrors, the 
SE rate may be significantly suppressed if the distance be- 
tween the plates is shortened below a critical value ]@]. Sup- 
pression of SE rate was observed by the first time in 1985 by 
Hulet et al UtIi by using a beam of excited Cesium atoms 
passing between mirrors and later on by other groups ifisl - 
I20I1 . There are excellent reviews on how strongly the envi- 
ronment can modify the SE rate of an emitter (such as atoms, 
molecules, or quantum dots), see, for instance, Refs lfl6l I21I - 

m. 

In the last decade the research on the radiative properties 
of atoms and molecules has received renewed interest due to 
progress in near-field optics. Indeed, advances in nano-optics 
have not only allowed the increase of the spectroscopical res- 
olution of molecules in complex environments ll26ll but have 
also led to the use of nanometric objects (e.g. nanoparticles 
and nano-tips) that to modify the lifetime and to enhance the 
fluorescence of single molecules ll27l l28ll . Another important 
example is the development of nanoantennas that can enhance 



the local optical field and modify the fluorescence of single- 
emitters ll29l[30ll . 

The advent of plasmonic devices has also opened new pos- 
sibilities for tailoring the emission properties of atoms and 
molecules. When a quantum emitter is located near a plas- 
monic structure it may experience an enhancement of the lo- 
cal field due the excitation of a plasmonic resonance, which 
affects the lifetime of an excited state. This effect has been 
exploited in the development of important applications in 
nanoplasmonics, such as surface-enhanced fluorescence and 
surface-enhanced Raman scattering ll3lij33ll and the modifi- 
cation of two-level atom resonance fluorescence ll34ll . It is 
also very important to point out the role that metamaterials 
may have on quantum emitters' radiative processes. Metama- 
terials are artificial structures with engineered electromagnetic 
response that may exhibit unusual properties such as negative 
refraction llssll . resolution of images beyond the diffraction 
limi t IBql, o ptical magnetism ]i;37..381, electromagnetic cloak- 
ing i39l koll and support of slow light propagation iiUl . Al- 
though the influence of some of these properties, such as nega- 
tive refraction, on the SE of atoms has been investigated ll43ll . 
many others remain unexplored so far. Recently, nanostruc- 
tured media with hyperbolic dispersion have emerged as a new 
class of metamaterials with many applications in biosensing, 
subwavelength imaging and waveguiding lliill . In addition, 
hyperbolic metamaterials exhibit a broadband enhancement 
in the electromagnetic density of states that has been demon- 
strated to have an important impact on the SE of molecules 
and quantum dots, allowing for several applications in quan- 
tum nanophotonics ll42ll44ll 

The aim of the present paper is to investigate the emission 
properties of an atom in the presence of a cloaking device. 
We specifically examine the behavior of the SE decay rate of 
an atom near a plasmonic cloak. The concept of plasmonic 
cloaking, theoretically proposed by Alii and Engheta ll46ll . is 
based on the scattering cancelation technique, in which a di- 
electric or conducting object can be effectively cloaked by 
covering it with a homogeneous and isotropic layer of plas- 
monic material with low-positive or negative electric permit- 
tivity. In these systems, the incident radiation induces a local 
polarization vector out-of -phase with respect to the local elec- 
tric field so that the in-phase contribution given by the scatter- 
ing object may be partially or totally canceled \A6, 47]. The 
first experimental realization of this idea was recently imple- 
mented for microwaves ll48ll . paving the way for many appli- 
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cations in camouflaging, low-noise measurements and non- 
invasive sensing. Here we focus on a novel application of 
plasmonic cloaking in atomic physics, demonstrating that the 
spontaneous emission rate of an atom can be drastically mod- 
ified in the presence of a plasmonic cloak. To the best of 
our knowledge we show for the first time that, in the dipole 
approximation, where light scattering from realistic particles 
is identically zero, the atomic SE rate reduces to its value in 
vacuum even for small distances between the atom and the 
cloak. This result not only proves that atomic SE decay rate is 
strongly modified by invisibility of the sphere but also sug- 
gests that the emission properties of an atom could be ex- 
ploited to probe the efficiency of a given plasmonic cloaking 
device. 

This paper is organized as follows. In Sec. HI] and Hill we 
describe the model and the calculation of the atomic SE rate 
whereas in Sec.|lV]we present, interpret and discuss our re- 
sults. Finally, Sec.|V]is devoted for our conclusions and final 
remarks. 



n. THE SPONTANEOUS EMISSION RATE OF A 
TWO-LEVEL SYSTEM 

In order to establish basic concepts and notation, as well as 
a convenient expression for the SE rate of an atom in the pres- 
ence of an arbitrary arrangement of bodies in its surroundings, 
we start discussing the atom-field dynamics. For simplicity, 
we consider an atom whose dynamics can be well described 
by two of its eigenstates. In the absence of interaction we as- 
sume that the lowest state 1 1 ) with energy Ei = — huiQ / 2 has a 
very long lifetime and a well-defined parity, while the highest- 
energy state 1 2) with energy E2 = fiujo/2 has opposite parity 
and a non- vanishing electric dipole coupling to Besides, 
the influence of the aforementioned bodies on the atomic ra- 
diative properties are taken into account in this model by the 
boundary conditions (BC) imposed by them on the quantum 
electromagnetic field modes. Bearing this in mind, let us 
assume that the system can be described by the well-known 
Hamiltonian iSotl 



H — Hax + T~L{ + Hint , 



(1) 



where T-Lat, 'Hf and Hm the atomic, the electromagnetic field 
and the interaction Hamiltonians, respectively. More specifi- 
cally, 

A, = ^hLJo{\2){2\ - |1)(1|) = ^huoa, (2) 

is the pure Hamiltonian describing the atomic internal dynam- 
ics, oja is the associated transition frequency and d-x,a-y,a-z 
are the Pauli operators. The second term in ^ is given by 



1 

8^ 



(3) 



where E(r) and B(r) are the quantum transverse electric and 
magnetic field operators. Furthermore, the atom is assumed 
to be much smaller than all other relevant length scales of 



the problem, allowing for the electric dipole approximation 
so that the third term in ([T]i reads 



-d-E, 



(4) 



where d is the atomic electric dipole operator. In the two- 
level approximation it takes the simple form d — d2iax, 
where d2i = (2|d|l) is the transition electric dipole moment 
that can be made real by an appropriate choice of the relative 
phase between the ground and the excited states. We also as- 
sume the atom has no permanent electric dipole moment so 
that (l|d|l) = (2|d|2) = 0. 

The electromagnetic and interaction terms may be recast 
in a more convenient form by expanding the electromagnetic 
field in their normal modes IdOII 

E{r)^iY,V^^C [acAc(r) - a^AJW], (5) 



Vx[acAc(r)-aU^(r)], (6) 



where the label C represents an arbitrary complete set of 
quantum numbers, ujq are the related eigenfrequencies, and 
ciQ , are the annihilation and creation operators that contain 
all quantum properties of the field and satisfy the commuta- 




tion relations [a, 



Note that A^(r) are classi- 



cal functions determined by the Helmholtz equation (V^ + 
A:|)A^(r) and the Coulomb gauge V • A^(r) = with 
the appropriate BCs. Substituting Eqs. ^ and ^ into (O 
and choosing an orthonormal set of functions Ai^{r), namely, 
J d'^r A*, (r) ■ A^;' (r) = (5^^/, one can see that up to zero-point 
contributions that do not contribute to the dynamics, we get 



Analogously, we have 

Hint = ^ift 

where we have defined 

9C - 



2'KUJr 



d2i-AJ(r) 



(7) 



(8) 



(9) 



In order to obtain to the SE rate, we have to determine the 
time evolution of the atomic energy average. Formally, this 
can be accomplished in the Heisenberg picture by writing an 
expression for a, (t) and then taking its mean value in the ini- 
tial state, where the atom is assumed to be in the highest en- 
ergy level and the field in the vacuum state. Applying the 
Heisenberg equation, we obtain 



C 



C 



,(10) 



(11) 



(12) 
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where we used the normal ordering for operators Sij and a|t. 
Also, (7, (jt are the Pauli lowering and raising operators for the 
atomic states and are related to by the following commuta- 
tion relations 

[<7,a,] =2<7, [<7^a,] = -2at, (13) 
[<7,at]=-a,. (14) 

The coupled Eqs. ( fTOl i. ( fTTI ). (fT2] i are quite difficult to solve, 
even for the simple system considered here. Besides, a direct 
integration in Eq. ( fT2l) shows that the relationship between 
ac^(t) and a^ify is non-local in time, with the field (atomic) 
operators depending on the atomic (field) operators in earlier 
times. Fortunately, considering that the atom is weakly cou- 
pled to the field, it is possible to simplify considerably these 
equations by performing a Markovian approximation and ar- 
rive at EEil 



where aj^(t) = a^(0)e 
lution and 



+ a(t)rK-^o)], (15) 
-%u)(^t jg j.jjg £j.gg ijgjj homogeneous so- 



V 



— \ — TTi6(x), 



(16) 



with V denoting the principal Cauchy value. Using the last 
two equations we can rewrite ( fTTT i in a simpler way 

C 

-2^[l + a,(t)]^|5cl'^(^C-^o)- (17) 

C 

Taking the average value of Eq. ( fTTI ) in the initial state |i) = 
2) (g) |vacuum) and using aj^(<) [vacuum) = 0, we obtain 

^r2i(r)t 



(CT.(i))2 = -l + [('j40))2 + l]e 
= -l + 2e-i^-M*, 



(18) 



where ((...))2 = (2|(...)|2) and we have used that {a-z{0))2 = 
1. The most relevant quantity in Eq. ( fTTI ) is, however, the 
spontaneous emission rate (r), given by 



r2i(r) 



Eld 

c 



21 



Ac(r)|2<5(c.^-c^o). (19) 



Equation (T9[ is the well known result for the SE rate of a two- 
level atom in the presence of bodies of arbitrary form. The 
influence of such bodies is coded in the functions (r) that, 
as mentioned before, are classical solutions of the Helmholtz 
equation. This is extremely convenient, as it makes the prob- 
lem amenable to various well-developed analytical and nu- 
merical techniques. Finally, as an easy but important check, 
and also for future convenience, one may use Eq. (T% to ob- 
tain the SE rate in free space. In this case, functions A.^ (r) 
are simply given by 



Ac(r) 



(r) = 



(27r)3/- 



(20) 



where Ckp is a polarization vector. Inserting Eq. ( l20b into Eq. 
( [19] ) we obtain 



r 



(o)_ 47r^a;o 



21 



„ 1 o J 



(21) 



3 h ' 

which is the well-known result firstly obtained by P.A.M. 
Dirac in 1927 JH. 



III. 



INFLUENCE OF A SPHERICAL PLASMONIC CLOAK 
IN THE SPONTANEOUS EMISSION RATE 



In this section we shall calculate the SE rate of an atom 
near a plasmonic cloak. Following the pioneering work of Alu 
and Engheta ll45[|46ll . the cloak is composed by a sphere with 
inner radius ai , permittivity ei (w) covered by a spherical shell 
with outer radius 02 > oi and permittivity 62 (w), as shown in 
Fig. [T] Both the core and the shell are nonmagnetic so that 




Figure 1. Cross section of a our spherical object showing the two 
radii ai and a2 > ai. The inner sphere and the covering layer are 
assumed to be made of isotropic homogenous materials, with permit- 
tivities and permeabilities ei (a;), /ii (a;) and £2(1^), ^2(1^), respec- 
tively. 

In order to obtain the field modes Akp(r) we consider the 
standard Mie approach for the scattering of a single plane 
wave (polarized in the e[ direction and propagating in the eg 
direction) by the plasmonic s phe re. The resulting field at an 
arbitrary observation point is i52ll53ll 



Akp(r')-A|,7)(r') + A(;^)(r'), 



(22) 



where A^^^'^\r') is given by Eq. ( |20] | and the scattered con- 
tribution can be written as a sum of TM*" and TE*" spherical 
waves 

1 



(27r)3/2 



v'xr(;^)(r') 



-V'xV'xG(;^)(r') 



(23) 



with 



kp 



(r') 



1=1 



i-i 



i2l + l)i' 
2l{l + l)u 



cf^kr'h\^\kr') 



(24) 
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OO 



(2/ + l)i' 



1=1 ^ ' 



cr'kr'h\'\kr') 



(25) 



In Eqs. (|24] | and dZSl l 6*' and Lp' are measured from the eg and 
e[ directions. yj™(6'',(^') are the spherical harmonics lf54ll . 

and ji{x) and /ip''(.T) are the spherical Bessel functions and 
the spherical Hankel functions of the first kind, respectively 
llsill . The coefficients and c[^^ are the so-called Mie co- 
efficients lf52l I53I1 . and here they are completely determined 
from the BCs at r = ai and r = 02. Also, the above expres- 
sions have a direct interpretation; the total scattered field is a 
superposition of the fields produced by all multipoles induced 
on the coated sphere. 




Figure 2. (Color online) Coordinate system used throughout the pa- 
per. The atom is placed at a distance r from the center of the plas- 
monic cloak. 

In order to calculate explicitly the sum in Eq. (fT9] l we now 
have to rewrite Eqs. (l24l) and (1251 ) in the unprimed coordinate 
system shown in the fig. |2] where the observation point (the 
position of the atom) lies within the eg direction. Formally, 
this can be accomplished by rotating the spherical harmonics 
in Eqs. (l24l i and (l25l l using the expressions 



'^"^"-dl^^,{h)Yr{0,v). (26) 



just have to evaluate the curls in (l24l i and (l25T l and reinterpret 
the relative angles 9' and 0' in terms of the spherical coordi- 
nates defined in the Fig. |2]as follows: r' ^ r, 9' ^ —9 and 
If' 0. Analogously, all the above equations and analysis re- 
main valid in the case in which the incident field is polarized 
in the e'2 direction and propagating in the eg direction; hence 
it suffices to change the relative angle ip' by tt/2 — <~p' . 

With those remarks in mind it can be shown that, by re- 
taining only the electric dipole contributions to the scattered 
fields, the electromagnetic field modes read 



.TM 



2(27r)3/2 



h\^\kr) 



+ hf\kr) 



'(r)= 



3ci 



.TM 



sin 26* 



2(27r)3/2 



hi'\kr) 



h^\krl 
kr 



''0] 



h[^\kr) 
kr 



h'[^\kr) cos^ 



(27) 



efc 



(28) 



where the prime in h'^\x) denotes differentiation with re- 
spect to the argument. 

It is worth emphasizing that we have two different approx- 
imations here: on the one hand, we have the (electric) dipole 
approximation for the atom, meaning that Aq = l-KcjujQ 3> 
ao, where ap is the Bohr radius; on the other hand, we have the 
(electric) dipole approximation for the cloaked sphere, which 
means that Aq = 2'kc/ijJq 3> a2- The former allows one to 
simplify the interaction part of the Hamiltonian, while the lat- 
ter justifies the neglect of all the scattering coefficients but 
^TA/ rpjjg (jipoie approximation for the coated sphere also im- 

pUesthat \c™\ < 1 and therefore |Aj;;"^(r)| < |A|^;"V)|. 



so that by substituting Eq. ( |22] | into Eq. (fT9] l we may keep 



only linear terms in a''^^\v), obtaining 



r2i(r) 



(29) 



kp 



where d\y^ ±i{S) are the so-called rotation matrices and au, f3k 
(and also 7^ = 0) are the respective Euler angles. However, 
in our case this is not necessary since we are rotating the coor- 
dinate system in Eqs. ( |23T l. (l24l i. and ( |25] ) by precisely the rel- 
ative angles 9' and ip' . Then, by symmetry, it is clear that one 



The first term in the previous equation gives the free space 
contribution (1211 1. while the second one is precisely the cor- 
rection due to the presence of the coated sphere. A direct 
substitution of Eqs. jZTl i and ( |28] | into Eq. (|29] | yields, after a 
lengthy but straightforward calculation. 



Ar2i(r) 

r(o) 
^ 21 



kor 



sin(fco?') 2 cos(fco?') sin(A:o7') 



{korf 



{kory 



knr 



2n'i{kor) 



sin{kQr) cos{kQr) sm{kQr) 



{korf 



(kor) 



kor 



(30) 



where Ar2i(r) = r2i(r) — ^21 ^ = Stt/Aq and ni{x) is the spherical Neumann function llsill . In order to obtain the 
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previous result we considered an isotropic atom, |d2ia:P = 
|c?2ij,P = M2i2p = |d2i P/3, and neglected losses by as- 
suming real permittivities at the transition frequency. 

Equation dSOb is the central result of this section and it 
demonstrates that SE rate of a two-level atom placed in the 
vicinity of a spherical plasmonic cloak is proportional to the 
first TM harmonic of the scattering coefficient of the coated 
sphere, |cf^^|. As shown by Alii and Engheta, this coeffi- 
cient, which corresponds to the electric dipole radiation, can 
be made to vanish by a judicious choice of material parame- 
ters ll45lJ47ll . As a consequence, the scattering cross section 
corresponding to the coated sphere will be greatly reduced 
largely dominates the scattering pattern, making 



smce c 



.TMl 



the system practically invisible to the incident radiation. This 
result has an interesting impact on the atomic SE rate, as it 
wiU be discussed in the next section. 



IV. DISCUSSIONS 

In the electric dipole approximation, the coefficient c[^^ 
reads 



TM 



.u- 



TM 



TM 



(31) 



where in the limit of low losses U^'^^ and are real func- 



rTM 



tions and can be written as 



TM 







1 1 -1 





[ka2 


? 







3 




7-1 -7' 


1 






2(7e2)'i 7'e2"' 2e^ 




11 -10 




1 




e,-i 




ka2 


-Y -1 






2(762)-^ l^£2^ £2^ 





,(32) 



rTM 



with 7 = ai /a2 and we made the vacuum electric permittivity 
equal to one, Eq = 1. 

The effectiveness of a cloak is related to the amount of scat- 
tering that is suppresed. In particular, an ideal cloak is char- 
acterized by no scattering at all, i.e., by the vanishing of 
It can be shown that a sufficient condition for that to happen 
isHlil, 



^ Ol ^ 3/ (£2-l)(2£2 + gl) 

^ a2 V (e2-£i)(2e2 + l)' 



(34) 



The complete set of material parameters for which invisibility 
can occur has been carefully discussed in Refs. ll45l - l47ll . and 
here we explore the case where j (ojo) < 1 < ^(ujo), with 
«,J = 1,2. 

In Figure[3]we show Ar2i (r) /T^2i ^ '^^e difference between 
the SE rate of a two-level atom near the coated sphere and its 
vacuum value normalized by Fj^', as a function of 02 /ai for 
a given distance between the atom and the sphere. We fixed 



l^Y^^ 0.006- 

r^™ 0.004- 

0.002- 
0.000- 
-0.002- 
-0.004- 
-0.006- 



= 1.51 
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Figure 3. (Color online) Relative spontaneous emission rate 
Ar2i(r)/r2i' as a function of 02/(11 for a given distance between 
the sphere and the atom. The sphere parameters have the following 
fixed values: ai = Aq/IO and £1 = 10. The blue dashed line and 
the red solid line show the results for £2 = 0.1 and £2 = 0.5, respec- 
tively. Vertical arrows indicate the values of 02 /ai for which the SE 
rate reduces to its value in the vacuum. For comparison, the green 
dotted-dashed line shows the result for the case where the coavering 
layer is made of the same material as the inner sphere (no cloak). 



the physical parameters characterizing the core by setting the 
electric permittivity of the inner sphere as £1 ~ 10 and, in or- 
der to be consistent with the dipole approximation, we chose 
the radius ai = Aq/IO. For comparison, in Figl3]also it is also 

plotted the result for AF21 (r) /T^2i '-^^^ where there is 

no cloak at all (green dot-dashed line), i.e., the relative SE rate 
in the presence of a single sphere with £1 = 10 as a function 
of its radius. In this case, as the radius of the sphere grows 
the SE rate deviates increasingly from its value in vacuum 
T^2i ' expected. In contrast, when the shell is taken into ac- 

j-^^) count, AF2i(r) approaches zero, and hence F2i(r) ^^21 ^ 
as the radius 02 of the cloak is increased. More interest- 
ingly, the SE rate F21 (r) is identical to its value in vacuum for 
02/01 ~ 1.09 (when £2 = 0.1) and for 02/01 ~ 1.51 (when 
£2 = 0.5), precisely the cases where condition ( |34] | is satis- 
fied for the respective choices of £2. As far as we know, this 
is the first situation where the SE rate of a quantum emitter is 
unaffected by the presence of a surrounding body for a broad 
range of separation distances. For values of 02/01 larger than 
the condition of invisibility, the contribution of the shell to the 
scattered field is dominant and the system becomes visible, so 
that the relative SE rate increases again as can be seen from 
Fig.S 

In order to investigate the dependence of the SE rate on the 
distance between the atom and the spherical plasmonic cloak, 
in Fig. |4]we exhibit a three-dimensional plot of Ar2i (r) /Fji^ 
as a function of both 02/01 and k^r. The parameters charac- 
terizing the inner sphere are the same as in Fig. [3] and the 
electric permittivity of the shell is £2 = 0.5. It is important to 
emphasize that the SE rate exhibits by and large an oscillatory 



behavior with the distance r except for 02/01 = -^38/11 ~ 
1.51 (highlighted by the red line in Fig. |4|, where the invisi- 
bility condition ( l34b is fulfilled and the SE rate is always equal 
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to for any value of kor. 




Figure 4. (Color online) Relative spontaneous emission rate as 
function of a^/ai and kor. Parameters of the inner sphere are 
ai = Ao/10 and ei — 10, while the covering layer has £2 ~ 0.5. 
The red solid line highlights the value 02/01 — ^38/11 ~ 1.51, 
for which the invisibility condition l l34t is fulfilled. We note also that 
0-2/0,1 — 1 corresponds to the SE rate of an atom in the presence of 
a dielectric sphere without cloak. 

To further explore the dependence of the SE rate upon the 
shell radius, in Fig. |5]we plot Ar2i (r) /Fj""* for different val- 
ues of the ratio 02/01; all other parameters being the same as 
those used in Fig. |4l Again, it is clearly seen that the pres- 
ence of the covering layer significantly reduces the amplitude 
of oscillation with respect to SE rate for a single sphere. The 
black dashed line corresponds to 02/01 = 1.40 and demon- 
strates that a reduction of about 80% for Ar2i (r) /T^^^ (com- 
pared to the case without the shell) can be obtained even 
with plasmonic cloaks for which the invisibility condition ( |34] | 
is not strictly satisfied. The green shaded area highlights 
the interval of possible values of a-ijai for which the rela- 
tive SE rate is reduced by at least 95% due to the inclusion 
of the shell. More specifically, the results show that when 
1.486 < 02/01 < 1.536, then F2i(r)/r^°^ is smaller than 5% 
of its bare value, regardless of the distance between the atom 
and the center of the cloak. This shows that a tolerance of 5% 
in the efficiency of the plasmonic cloak is possible in a rela- 
tively wide range of ratios 02/01, highlighting a serendipitous 
robustness of our results. 

Let us now consider some realistic parameters in order to 
assess the experimental viability to test the results here dis- 
cussed. A well-suited quantum emitter is a Rubidium atom 
prepared in a Rydberg state with principal number n = 51 
and magnetic number m = 50, since it is possible to show that 
such an atom can be well described by a two-level system and 
therefore our previous discussion can be applied pSf). More- 
over, the corresponding transition frequency is vq = 2ttujq « 
51.099 GHz, well in the range of recent experiments with 
plasmonic cloaking devices ||48[ l56ll . Suppose the Rubidium 
atom is placed near a spherical cloak composed of a non- 
magnetic sphere of radius oi and permittivity ei = 10 cov- 
ered by a spherical shell with outer radius 02 and permittivity 
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Figure 5. (Color online) Relative SE rate as a function of kor for 
different values of a^jox. The blue solid line correspponds to the 
absence of a cloaking shell, namely, 02/0-1 = 1. The black dashed 
line corresponds to 02/01 — 1.40 and shows that the covering layer 
significantly reduces the amplitude of oscillation. The red dotted- 
dashed line shows the case where the condition i34i is fullfilled, 
02/01 = ^38/11 ~ 1.51. The shaded area corresponds to the 
interval 1.486 < 02/ a\ < 1.536 and highlights the range of values 
of 02/ai for which Ar2i(r)/r2i' is reduced to 5% or less than its 
value in the presence of a single sphere. 

£2 = 0.5, being both permittivities evaluated at the transition 
frequency ujq. As before, let us take oi = Aq/IO = 587.1/im. 
Therefore, the ideal value for 02 for a perfect cloaking is 
02 = ^38/11 X 587.1 ^ 887.5//m. If, instead of per- 
fect cloaking, we require only that the SE rate is suppressed 
by at least 95%, the possible values of 02 are in the interval 
872.4/xm < 02 < 901.8/iTO. In other words, one could vary 
02 by ztlbfim around the ideal value and would still have 
a very efficient reduction of the SE rate (> 95%). Had we 
been interested in a 99% of the SE rate, the allowed range for 
02 would be narrower, namely, 884.8/im < 02 < 890.7/im 
(in this case, 02 could be changed by ±3/im around the ideal 
value). These numbers show that, at least in principle, our re- 
sult could be tested with current apparatuses and techniques. 

V. CONCLUSIONS 

We have investigated the SE rate of a two-level atom placed 
in the vicinity of a plasmonic cloak composed of a coated 
sphere. In the dipole approximation, i.e., when the coated 
sphere is much smaller than the atomic transition wavelength, 
we have demonstrated that the difference between the SE rate 
in the presence of the cloak and its vacuum value is directly 
proportional to the first TM scattering coefficient. Since this 
coefficient can vanish under certain circumstances, as estab- 
lished by Alu and Engheta [46], this result implies that the SE 
rate of a two-level atom can be identically reduced to its vac- 
uum value. We have analysed the dependence of the SE rate 
on the distance between the atom and the plasmonic cloak; we 
conclude that the reduction of the SE rate is significative for a 
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large range of distances, even for small distances between the 
atom and the cloak. We have also investigated the dependence 
of the SE rate on the geometrical parameters of the cloak, such 
as the ratio between its inner and outer radii, as well as on its 
material parameters, such as the electrical permittivities; we 
find that the strong suppression of the SE rate is robust against 
the variation of both the geometrical and material parameters 
of the cloak, as a manifestation of the robustness of the plas- 
monic cloaking mechanism itself. 

In order to envisage the possible experimental verification 
of our findings, we have made realistic estimates for the SE 
rate of atoms in the vicinities of a plasmonic cloak. Since 
the frequency operation range of state-of-the-art plasmonic 
cloaks is typically of a few GHz ll48l l56ll . we argue that the 



presented results could be experimentally verified as there are 
many atomic species that spontaneously emit in this frequency 
range, such as Rubidium, the explicit example we have con- 
sidered. As a result, we suggest that the observation of the 
reduction of the atomic SE rate to its vacuum value in the pres- 
ence of plasmonic cloaks could be explored as an alternative, 
quantum probe of the effectiveness of these devices. 
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